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Abstract 



Extraction of spin observables from vector meson photoproduction on a nu- 
cleon target is described. Starting from density matrix elements in the vector 
meson's rest frame, we transform to spin observables in the photon-nucleon 
cm. frame. Several constraints on the transformed density matrix and on 
the spin observables follow from requiring that the angular distribution and 
the density matrix be positive definite. A set of constraints that are required 
in order to extract meaningful spin observables from forthcoming data are 
enunciated. 
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I. INTRODUCTION 



In an earlier paper we emphasized that the angular distribution of the pseudoscalar 
mesons that arise from the decay of photoproduced vector mesons does not depend on the 
vector meson's vector polarization, but only on its tensor polarization and that standard 
single and double spin observables need to be defined in the overall photon-nucleon center of 
mass frame. It was also found that a simple description for the decay angular distribution 
in the 7iV cm. frame is obtained by using the angle between the decay meson's velocity 
difference vector and the direction of the photoproduced vector meson. The main purpose 
of this paper is to formulate a procedure for extracting meaningful spin observables from 
the analysis of forthcoming vector meson photoproduction data to allow one to examine 
conventional spin observables. These spin observables are subject to known rules concerning 
their forward and backward angular behavior p[. The nodal structure of spin observables, 
e. g., their production angle dependence, might reveal important underlying dynamics such 
as baryonic and mesonic resonances. 

Here we show how to extract meaningful spin observables under the assumption that 
analysis of the photoproduction of vector mesons will yield a vector meson rest-frame den- 
sity matrix. Since new data are not yet available, we invoked older 1968 Aachen et al. 
information |^ and found that some of their vector meson rest-frame density matrix results, 
when transformed to spin observables in the photon-nucleon center of mass frame, violated 
basic constraints and therefore need to be rejected. The grounds for that rejection was that 
some of their elements, even including their stated uncertainties, yielded non-positive, and 
therefore unacceptable angular distribution functions. That observation, which we subse- 
quently found to be related to a set of constraints deduced by others earlier @J^, led us to 
examine the various constraints based on the positivity of the density matrix. 

In Sec. II, we analyze the limits on the tensor polarization provided by the simple re- 
quirement that the angular distribution of the decay mesons be a positive definite function. 
Simple limits on the tensor polarizations follow from evaluating the decay angular distribu- 
tion at selected angles. In Sec. Ill, the limits on observables due to positivity of the density 
matrix are discussed. In Sec. IV, constraints on the vector meson's density matrix due to 
Daboul's analysis using Schwarz inequalities, are invoked and analyzed. The Schwarz 
inequalities described by Daboul can be re-expressed as four separate conditions on spin 
observables. However, all but two of these conditions are already contained in the simple 
requirement that the angular distribution should be positive definite. The two remaining 
conditions involve not only the tensor polarization, but also the vector meson's vector polar- 
ization. These conditions could also be used to limit the vector meson's vector polarization. 
Finally, in Sec. V the method for extracting spin observables from actual data is outlined. 

Use of these basic constraints should be included in the fitting procedure to assure that 
general requirements concerning the angular dependence of spin observables, especially at 
forward and back angles, are satisfied and could then be used to deduce interesting new 
dynamics. 
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II. THE DECAY DISTRIBUTION 



In photoproduction of a vector meson (p, (p) on a nucleon target the final vector me- 
son decays into two pseudoscalar mesons. The angular distribution of the decay provides 
information about the spin-state of the vector meson. However, only information about 
the tensor polarizations T^,T2i,T^ can be obtained. The angular distribution of the pseu- 
doscalar decay mesons is given by 



^^vie) [1-V2 TUe, $) cum 



= -^^viO) [1 - ^-T^o{3cos^e-1) + V3T^i sin 2^ cos 0- 73 sin2^cos20]. (i; 

Here = y^^2/i is a spherical harmonic function and the angles 6, refer to the direction 
between the velocity vector difference, Av = vi — V2, and the momentum vector of the 
produced vector meson; Vi and V2 refer to the velocity vectors of the two decay mesons in 
the overall photon-nucleon center of mass frame. Use of these angles simplifies the expression 
for the angular distribution in the overall photon-nucleon center of mass frame in which spin 
observables are defined. Note that the spin observables depend on the vector meson 
production angles 0, $, as well as on the total cm. energy. The factor ^, which arises from 
describing the vector meson decay in the overall center of mass system and from a density 
of state factor, is given by: 

(sm u + (^) cos^ a)^!'^ 

where mp, Ep are the vector meson's mass and energy. 

The decay angular distribution {6,(j)) does not depend on the vector meson's vector 
polarization and as shown above includes only the vector meson's tensor polarization. Once 
the angular distribution is measured and vector meson rest frame density matrices are pro- 
vided, it is necessary to map that data over to the angles 6, 0. One can then project out the 
vector meson's tensor polarization from the normalized ratio W^^j.^{9, (j)) = W {0, 0)/^, as 
described in Sec. V. 

The tensor polarization must take on values that allow the angular distribution function 
Vr^(i9, 0) to be positive definite. By selecting the angles ^,0 one can use that obvious 
condition to extract allowed limits for the tensor polarization.^ In the first and second 
columns of Table I a list is given of specific choices of angles 9, and the resulting conditions 
on T2Q,T2i,T^. Such constraints also arise from direct conditions on the density matrix, as 
will be seen in Sections HI and IV. 

Thus from Table I, we see that the simple requirement that the decay angular cross 
section be positive yields limits on the possible tensor polarization. We now consider other 



The allowed ranges for the tensor polarization can also be deduced by considering the spin-state 
occupation amplitudes in the pure state limit. 
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ways to recognize constraints on the spin observables and the associated density matrix. In 
the next section, we describe the constraints on T^,T2i,T^ that follow from the positivity 
of the density matrix. 



III. LIMITS ON OBSERVABLES FOR A POSITIVE DEFINITE DENSITY 

MATRIX 

Recall that for a general observable Q the classical ensemble average is 

<n>= ^^^^<^M^>^ (3) 

where Ua is the positive definite probability for finding a beam particle pointing in the 
direction stipulated by the Euler angle label a. Note that the above is a classical average, with 
the quantum effects isolated into the expectation value for each beam particle < a;|f2| a > . 
The spin density matrix of the vector meson is defined as 

P = \ > < o: \ , (4) 

a 

where Ua is non-negative. The helicity matrix elements of the density matrix are 

pxx' = ^ < A I a > ciJq, < a I A' > . (5) 

a 

The classical ensemble average for observable Q is now obtained from the density matrix p 
as 

The density matrix p is positive definite, which can be shown as follows. Let us define Q 
a set of vectors vx by its elements in a space 

=< a\X > ^JuJ^ . (7) 



The elements pw of the density matrix can now be written as dot-products in a space of 
the set of vectors vx , 

PAA' = (^^A,^^A')- (8) 

Now for any vector X 

X^pX = J2 XlpxX'Xy = (Z XxVx, E ^A'^^A') > 0, (9) 
A,A' A A' 

which displays the positive definiteness of p. 

At this point we explore the linear constraints on the matrix elements of p implied by 
Eq. (1^). The density matrix p is a 3 x 3 matrix with elements px,x' where the helicity A takes 
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the values 1, 0, —1. Since the production of the vector meson occurs via a parity conserving 
mechanism, the spin density matrix elements satisfy the symmetries 

Pw = P*\'\, and paa' = (-l)^"^'p-A-A'- (10) 
The density matrix takes the form 

/ (l-Poo)/2 ^piQ + i^pw pi-i \ 

P = ^Pio - i^Pio Poo -^Pio + '^^^Pio • (11) 

\ pi_i -^pio-i^pio (l-poo)/2 / 

In the language of spin observables Py , T^q, T^i, the density matrix can be written as 

p=i[/ + ^S^-P^ + r-Tn, (12) 

where S is the spin-1 operator and r is the symmetric traceless rank-2 operator with cartesian 
components Ty = ^{SiSj + SjSi) — 26ij. In matrix form this becomes 



1 

^=3 



^ ' y 2 20 2\l 2^ y' V 2 21 V'J-f22 



V2 2I 2V2V ''-'i//^V2 21 



(13) 



If the vector X in Eq. (|^) is such that the combination Xl^Xy is symmetric under the 
exchange of A and A', the obtained constraints are similar to the constraints derived in 
the previous section for a positive decay angular distribution. The constraints in that case 
involve only the tensor polarizations Tg^^jTj^, and Tg^, and not the vector polarization . 
This is because symmetric combinations X^Xa' only pick out the symmetric part of the 
density matrix paa', i- e., in this case the part that has even rank. The antisymmetric rank- 
one part of p, which is due to the vector polarization, then gives no contribution to X^pX. 
This is exactly the symmetry selection made if one considers the decay angular distribution 
of Section II, and which was discussed in detail in Ref . |I| . The resultant linear constraints 
are listed in columns 1 and 3 of Table I. 

Even relations involving the vector polarization Py can be obtained from Eq. (|^), if 
the above symmetry restrictions are not invoked on X^^Xa'. However, the resulting linear 
constraints involving P^ are only of academic interest because P^ cannot be measured from 
the decay angular distribution. These additional constraints are therefore not listed in this 
paper. 



IV. SPIN OBSERVABLE LIMITS FROM SCHWARZ INEQUALITIES 

Additional constraints on the density matrix are obtained using Schwarz inequalities as 
described in Ref. Q. Namely, from the previously derived Eq. (^ and 

\{vx,vxi)\ < \vx\\vxi\ , (14) 

follows 
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Ipaa'I < a/Paa Pa'a' • (15) 

Similar constraints exist for differences or sums of matrix elements p\\i. Such constraints 
were exploited in Ref. [Q. In this case, two additional inequalities can be derived 

Ipaa' + P-AA'I < y^2(pAA + Pa-a) PA'A', (16) 

and 

Ipaa' - P-AA'I < y^2(pAA - Pa-a) PA'A'- (17) 

From Eqs. (|l5| , p!6| , |T7D one finds several quadratic constraints. Using Eq. (|10D and the 
property that the diagonal matrix elements Pi,i,Po,o; and P-i-i are non-negative or 

1 + ^ Tj^ > 0, (18) 

1 - 72 T^^ > 0, (19) 
some of the Schwarz inequalities collapse to linear conditions 



1 + ^1^20 + > 0, 



(20) 



1 + sI\tI - VST^, > 0. (21) 
However, one also finds two very useful quadratic conditions that involve the squares of 

9{P^f + 12{T^,f < 8(1 - V2T^,){1 + ^T,^), (22) 



and Tj^, 



9{P^r + 12(T,\)^ < 4(1 - V2T^,){1 + ^T^, - VST^, 



(23) 



The resulting restrictions on the observables are listed in columns 1 and 4 of Table I and 
Table II. As one can see, the obtained linear rules are equivalent to those discussed earlier. 
Again we omit from Table I linear conditions that include Pj^. 
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V. DATA ANALYSIS METHOD 



The Aachen et al. 0] collaboration measured the reaction 7 + p —>■ + p using an 
unpolarized photon beam and an unpolarized proton target. The final proton recoils when 
a p meson is produced and the p subsequently decays into 7r+ and vr" mesons, both of which 
are detected. Hence, the angular pion distribution in the final state is measured. A fit 
to this angular distribution in the p meson rest-frame yields three density matrix elements 
from which we can reconstruct their pion angular distribution in the p meson rest frame. 
This reconstructed decay pion angular distribution is called W^{9, 0), where 9, define the 
direction of Tr"*" in the p rest frame. 

The above pion angular distribution depends on the spin state of the produced p meson, 
which is described by a 3 x 3 spin density matrix Pa,a'- The angular distribution only depends 
on the three real elements poo, ^Pio, Pi-i- Values of these elements in the rest frame of the 
p meson are published by Aachen et al. for a set of several photon beam energies and 
vector meson production angles. In order to study reaction mechanisms, however, we are 
interested in spin correlations (single and double spin observables) that are defined in the 
overall center of mass frame. How does one obtain spin correlations from these previously 
published density matrix elements? 

Our aim is to use the Aachen et al. data in the form of W^{6,(f)) to obtain the pion 
angular distribution in the photon-nucleon center of mass frame and re-analyze it in terms 
of the spin correlations. If one or more particles in the reaction are polarized [|I|, such 
future data can be analyzed in a similar way to extract meaningful spin correlations. Our 
procedure is therefore preparation for analysis of future experimental results from Thomas 
Jefferson Laboratory with polarized photons [0]. 

The first step is to obtain the angular distribution W^{6,(j)) in the p meson rest frame 
from the values of the elements poo? ^Pw, Pi-i using 

W^{e,^) = —[ ^^ + ~ ^ cos^ 9-V2 ^pio sin 29 cos - pi-i sin^ 9 cos 20] . (24) 
47r 2 2 

Then one constructs the angular distribution {9, 0) in the 7-nucleon center of mass using 

(sm 9 + (^) cos^ 9y^''^ 

In the cm. frame the variables are 9 and 0, where 9, are the angles between the relative 
velocity of the two decay pions in the 7-nucleon cm. frame. These angles are related to the 
angles of the decay meson in the vector meson's rest frame by: 

E 

^ = arctan( — ^tan^), (26) 
Trip 

and (f) = (j). 

The next step involves including the known kinematic factor C,{9) in W {9, 0) see Eq. (|j). 
Aside from the overall factor of ^{9), the remaining angular behavior is expressed as a series 
in spherical harmonics: 
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W^{9,(l)) = ^^{9)[{1- JIt^q{3cos^9-1) + Vs TJ^ sin 2^ cos - sin^ ^ cos 20] . (27) 

47r V 2 

Next we define 0) using Eq. (g^) and W"" {9, 0) = e(^) x Wl,,^{9, 0). 

Finally, we project out the three spin observables, i. e., the three tensor polarizations 
of the vector meson, ^205^21' ^22' from w'^^^j^{9 , (p) using spherical harmonics Yim(^, 0)'s. 
These spin observables can then be studied as function of photon energy and the vector 
meson production angles O, $. 

Once the spin observables are properly defined and have correct production angle de- 
pendence, one can visualize the role of the tensor polarization in 3D displays of the decay 
angular distribution. Examples are given in Fig. 4 for the case of = 0°, 180° and pos- 
itive T20 and in Fig. 5 for the case of = 0°, 180° and negative T20; with in both cases 
T21 = T22 = 0. In Fig. 6, a more realistic case, based on Ref. for = 3GeV and = 70° 
is shown; namely, for T20 = — .72,T2i = — .21,T22 = .19. One can therefore associate a shape 
of W with each point in the allowed T2o,T2i,T22 space. 

In the process of carrying out these steps, we found four cases of the Aachen data that 
do not satisfy the constraints in Tables I and II. Therefore those sets had to be rejected. 
Subsequently we found another author had also rejected some of the data using similar 
general constraints. 

It would therefore be best to incorporate these constraints directly into the data analysis. 

VI. CONCLUSION 

For analysis of experimental data for vector meson photoproduction one should describe 
the resulting angular decay distribution in the overall cm. frame instead of in the 
commonly used vector-meson rest frame. In the cm. frame, one should use the angles 9, (j) 
of the relative velocity vector Av = Vi — V2- In the transformation from the vector- meson 
rest frame to the cm. frame, and due to the use of the angles 9, 0, a kinematical 
factor ^ of Eq. (|D needs to be included. Furthermore, constraints should be satisfied by the 
observables T2q,T2i,T^. All of these constraints can be derived from the positivity of the 
density matrix. 

We have also explored restrictions which follow from positivity of the eigenvalues of the 
spin density matrix p such as the conditions that det p > and Trace [p^] < 1, see Ref. 0. 
The relations obtained from these restrictions are respectively cubic and quadratic in the 
spin observables and involve the vector meson's vector polarization . A more complete 
set of relations can be obtained by exploring the explicit forms of the roots Xi,X2, of the 
eigenvalue equation of p. We impose the conditions that < Xi < 1, Xi + X2 + X3 = 1, 
and xf + xl + xl < 1. For example, one root is Xi = |(1 + :^^2^ + v^^2^)- From < 
Xi < 1 relations follow that are similar to the rules in Table I. The other two roots are 
^2 = 1(1- i73^2^o - fT^2) - \^Py + -zTl + {^T2,-^T22Y and x^ = l-x^- x^. Both 
X2 and involve Py and the above root rules lead to quadratic constraints that are included 
in Table II. 

The linear and quadratic constraints on Tg^j, Tg^, T22 mean that the allowed domain in the 
3-dimensional — T21 — space is confined. As examples Figs. (1-3) show the allowed 
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areas for the 2- dimensional subspaces — Tj^ , — , and — T^i respectively. 
The dashed lines indicate the various upper bounds and the solid lines represent the lower 
bounds associated with the various constraints. The allowed regions are shaded. Figs. (1-3) 
refer in each case to 2-dimensional subspaces and only those constraints are shown that are 
independent of the values of the observable that would play the role of the third dimension. 
The full 3-dimensional representation of constraints is richer. 

We close with some remarks about double spin observables. For a polarized photon beam 
again the angular distribution of the decay mesons can be measured. This decay distribu- 
tion now depends on the single spin observables ^20 '^21 '^22 well as on the double spin 
correlations C22o,C22i,C222,C^2i,C^22,C]2i,C]22 (Ref. 0). (Again the vector polarization 
and correlations with the vector polarization cannot be measured.) Similar to the method 
described in Sections II and III and based on positive decay angular distributions and posi- 
tivity of the complete density matrix for this case, linear and quadratic relations involving 
double spin observables can be derived. 

The constraints on the spin observables presented in this paper should be incorporated 
directly into the analysis of forthcoming data. 
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TABLES 



Constraints 


> 


X^pX > 


Schwarz 


1 - V2r2o > 


o = fA = l 


Poo > 

Pii >o 


Poo > 

Pii >o 


i + ^r2i>o 


e = h;^ = l 


1 - 2V23fipio > 
1 + 2V23fipio > 




l + ^T22>0 




1 + 2pi_i > 




1 + ^^20 + V3r22 > 

1 + ^T20 - V3r22 > 

2 - ^T2o - x/3T22 > 


e = fA = ^ 
e = f,^ = o 

- 4'^ - 2 


1 - Poo + 2/9i_i > 
1 - Poo - 2/91-1 > 


|pi-iP < piiP-i-i 
|pi-iP < piiP-i-i 
see caption(*) 


1 - 373^20 + 4^21 > 
1 - 373^20 - 71^21 > 


;^ = 1 


1 + Poo - 43fi/9io > 
1 + Poo + 43f?pio > 




i + W^2i-A^22>o 
i-W^2i-^r22>o 


^ = ^c; ^ = 
e = ec;^ = 7r 


1 - 43fipio - 2pi-i > 
1 + 43fipio - 2pi-i > 





TABLE L Linear Constraints for T2o,T2i,r22; 6c = arccos(;^). (*) This condition follows 
from Xi < 1. 
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Constraints 


Schwarz Inequality 


Other 


9P2 + UTl < 8(1 - V2T2o){l + ^T2o) 


IpioP < PiiPoo 




9P2 + 12r|i < 4(1 - V2T2o)(l + -^T2o - VST22) 


1/910 - P-lOp < 2|pil - Pl-l|/900 


X2 > 


9P2 + i2Tl < 8(1 + ^r2o)(2 - ^T2o + V3r22) 




X3 < 1 


9P2 + i2Tl + I2TI2 + 6r|o < 12 




rr[p2] < 1 



TABLE II. Quadratic constraints for Py,T2o,T2i,T22 



12 



FIGURES 



1 



T21 









0.5 








0.5 i 




-0 . 5 









4- T20 



FIG. 1. Allowed domain in the space spanned by T^q and Tj^ independent of the value of T22- 
Dashed lines are upper bounds. Solid lines are lower bounds. Linear constraints are from Table 1 
and quadratic constraints are from Table II. The shaded area represents the allowed region. 
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FIG. 2. Allowed domain in the space spanned by and T22 independent of the value of T^i- 
Dashed lines are upper bounds. Solid lines are lower bounds. Linear constraints are from Table 1 
and quadratic constraints are from Table II. The shaded area represents the allowed region. 
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FIG. 3. Allowed domain in the space spanned by T22 and T21 independent of the value of TgQ. 
Dashed lines are upper bounds. Solid lines are lower bounds. Linear constraints are from Table 1 
and quadratic constraints are from Table II. The shaded area represents the allowed region. 
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FIG. 4. Forward/Back W for T20 > 0. This is the distribution in terms of the angles 9, (/) with 
the up direction being the direction of the produced vector meson's momentum. 




FIG. 5. Forward/Back W for T20 < 0. This is the distribution in terms of the angles 9, (j) with 
the up direction being the direction of the produced vector meson's momentum. 
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FIG. 6. Shape of W for Tao = -.72, T21 = -.21, T22 = -19. This is the distribution in terms of 
the angles 0, (j) with the up direction being the direction of the produced vector meson's momentum. 
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